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Abstract — The everlasting bandwidth limitations in wireless 
communication networks has directed the researchers' thrust 
toward analyzing the prospect of overloaded Code Division 
Multiple Access (CDMA). In this paper, we have proposed a 
Genetic Algorithm in search of optimum signature matrices 
for binary-input synchronous CDMA. The main measure of 
optimality considered in this paper, is the per-user channel 
capacity of the overall multiple access system. Our resulting 
matrices differ from the renowned Welch Bound Equality (WBE) 
codes, regarding the fact that our attention is specifically aimed 
at binary, rather than Gaussian, input distributions. Since design 
based on channel capacity is computationally expensive, we have 
focused on introducing a set of alternative criteria that not only 
speed up the matrix formation procedure, but also maintain 
optimality. The Bit Error Rate (BER) and Constellation measures 
are our main criteria propositions. Simulation results also verify 
our analytical justifications. 

I. Introduction 

Since the dawn of synchronous Code Division Multiple 
Access (CDMA), there has been an everincreasing demand to 
improve their superior capabilities, such as channel capacity. 
Quality of Service (QoS) and the per-user data bit rate. One 
main obstacle in the path of system engineers, however, is to 
increase the number of CDMA channel users while sustaining 
reliable transmission. Of course, if the number of users is less 
than or equal to the available spreading factor, the utilization 
of completely orthogonal codes, e.g. Hadamard codes, leads 
to satisfaction in the matter. Nevertheless, no such orthogonal 
codes are attainable for an overloaded CDMA system; thus 
errors caused due to Multi Access Interference (MAI) as well 
as inherent system noise are inevitable. 

The analysis of the overloaded CDMA concept comes to life 
when channel bandwidth limitations are realized, therefore a 
great deal of research has been aimed at the subject in recent 
years. For example [l]-[2] have focused on pseudonoise (PN) 
spreading, [3] has presented PN/OCDMA (PN/0) signature 
sets and [4] -[6] have presented well-behaved signature sets 
for binary-inpuQ CDMA channels, namely COW matrices 
and [7] -[8] have evaluated the overloaded noiseless and noisy 
channel capacity for different orders of dimension. 

The renowned Welch Bound Equality (WBE) codes have 
been introduced in [9] and justified to reach maximum capacity 
for Gaussian input distributions [10]-[11]. However, regarding 

'By binary-input, we mean tliat the input vector to tlie system is assumed 
to be binary, i.e. itl's. 



the fact that this paper's attention is specifically on binary 
input signals, the capacity-reaching justifications for WBE's 
are no longer true. As a matter of fact, according to simulation 
results, WBE's exhibit diverse performances for binary inputs 
and in most cases below optimum characteristics. 

Among the literature, [4] -[6] are the only that deal with 
binary-input overloaded CDMA. These papers have proposed 
the construction algorithms for matrices with binary (—1 and 
1), ternary (—1, and 1) and generally limited elements in the 
CDMA signature sequences. Although limiting the elements of 
the signature matrix causes design simplicity, it usually leads 
to optimality limitations. Conversely, in this paper, we have 
not confined the elements of the optimum signature matrices, 
apart from the fact that they should be real. 

It is noteworthy that none of the binary-input signatures 
schemes that have been proposed in the literature (including 
the WBE and COW matrices) guarantee optimality regarding 
to channel capacity maximization in different SNPfl 

In this paper, we basically pursued the following objectives: 

• Evaluating the key equations relating any signature se- 
quence matrix with its corresponding binary-input syn- 
chronous CDMA channel capacity. 

• Introducing a set of appropriate alternative criteria, in- 
stead of the channel capacity itself, that could be opti- 
mized over signature matrices in order to simplify and 
speed up the procedure of attaining optimum capacity 
code sets. 

« Optimizing the proposed criteria using a real-valued Ge- 
netic Algorithm and comparing the results with analytic 
evaluations. 

In section we have presented the basic preliminaries of 
our research. The general definition of a binary-input AWGN 
CDMA channel, followed by an outline of our main concern, 
i.e., the channel capacity maximization, is presented. Hence, 
with the results of the formulations presented in this section, 
the precise channel capacity could be obtained for any specific 
signature matrix. 

In section Uni some criteria for the optimality of a signature 
code matrix are presented. These measures are proved to be 
in coherence with optimum (or in some simplifications sub- 
optimum) capacity characteristics and prepare the grounds for 

^Most of the work done in the Hterature assume noiseless and high SNR 
conditions rather than also evaluating their propositions in relatively low 
SNRs. 



rapid selection of the signature matrix set for the overloaded 
system. The main setback in choosing the perfect measure is 
that a compromise must be solved between each criterion's 
complexity and its optimality performance. Nevertheless, our 
proposed criteria not only maintain near-optimum characteris- 
tics, but decrease the computation order to a great extent. 

There are two main approaches that have been proposed in 
maximizing the overloaded sum channel capacity, the direct 
channel capacity maximization and the indirect Bit Error Rate 
(BER) minimization schemes. The indirect scheme is based 
on our conjecture that matrices that exhibit low BER char- 
acteristics, also behave near-optimally regarding to channel 
capacity maximization. Based on the second set of criteria, in 
subsection IIII-CI we have justified the proposition of a set of 
Constellation measure^ 

In section |IV] we have carried out the Genetic optimization 
Algorithm in order to execute and deeply compare all the pri- 
orly explained criteria optimizations. Finally the conclusions 
and future works are depicted in section |Vl 

II. Problem Definition and Preliminaries 

Synchronous CDMA systems, regardless of the fact that 
they are overloaded or not, are formulated in a consistent 
manner as follows. 



Y = AX + N 



(1) 



where A is an to x 71 matrix with CDMA signature columns. 
Due to the fact that this paper's concentration is basically on 
overloaded synchronous CDMA systems, the number of users 
(rt) is larger than the signature spreading factor (m). Also, 
X is the user column vector with entries {±1} and N is the 
AWGN vector N = [Ni,.. . ,iV,„]^, such that N^s are i.i.d 
random variables. 

We shall also denote the overloading factor of such a system 

as, 

/3=^ (2) 
m 

In general, the elements of A are free to take any real 
or complex value, but in our focus, A should be composed 
of real-valued elements with the only restriction of having 
normalized signature columns to promote fairness between all 
users. 

It is clear that the sum channel capacity for a given A matrix 
is as below, 

C(A) = maxI(X;Y) (3) 

p(X) 

Note here, that the maximization is over all possible mul- 
tiphcative distributions on X. However, according to the 
conjectures mentioned in [4] and [7], for every matrix A, the 
capacity is met for a uniform input distribution. Thus we have, 

C{A)^I{X;Y) , p{X^X)^l^ (4) 

^Constellation measures deal with the possible output projection points in 
the image field for noiseless conditions. 



In order to facilitate our path towards generating the optimum 
overloaded CDMA matrices, we need to initially evaluate the 
maximum channel capacity. We define the maximum channel 
capacity as. 



max C(A) 



C(n, TO, (Tat) 

Taking into account that A is deterministic, we have, 

I{X;Y) ^ h{Y) - h{N) 



(5) 



(6) 



Also, keeping in mind that the Gaussian noise elements are 
independent, fN{N) is as follows. 



fN{N) 



( 1 
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(7) 



Based on [7] and due to the uniform distribution of X, the 
probability distribution of Y can be computed as below. 
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(8) 

where and A^ are the i*'' entry of Y and the i*'' row of 
matrix A, respectively. In addition, the summation is over all 
possible X e {±l}"^i. 

According to dHJ, we can have 



hviY) 



/y(r)log2(/y(r))dyi...dy„ 



(9) 



Note that in order to evaluate C(A), we need to subtract 
h{N) from (|9]), where h{N) can also be computed using the 
definition of differential entropies. 

It's noteworthy to mention that using the relations men- 
tioned above, one can evaluate the sum capacity of any arbi- 
trary matrix, which is the basis of our comparisons throughout 
the paper. According to this achievement, we have the ability 
to compare various proposed matrices regarding to the opti- 
mality of their channel capacity. 

In the next section, we have proposed the different possible 
criteria and methods in choosing optimum signature matrices 
for different dimensions and SNRs. The schemes are presented 
and deeply evaluated under different scenarios, such as the 
Capacity, BER and Constellation Criteria. 

III. Optimum Matrices Criteria 

As priorly explained, we have tried to match the optimality 
of a signature matrix capacity with certain criteria, so that it 
would be possible to find an optimum matrix with arbitrary 
dimensions and SNR. This approach helps the system designer 
to easily implement an overloaded CDMA structure without 
the fuss of probing blindly into any feasible matrix. To put 
it straight, in this approach the system designer just has to 
optimize the certain criterion, while being sure of the sub- 
optimality (if not optimality) of the generated signature matrix. 



0.950 













— ►—Optimum for 5dB 






— 1 — Optimum for 8dB 






— Optimum for lldB 






— e — Optimum for Minimum BER 






The Extreme UpperBound 











'3 4 5 6 7 8 9 10 11 12 13 



SNK (clB) 

Fig. 1 . The Per-User Channel Capacity of Optimum Code Matrices (Capacity 
and BER) with ^ = | Versus Different SNRs 

A. The Capacity Scenario 

Now that we have provided the grounds for evaluating the 
sum channel capacity of an overloaded synchronous CDMA 
system, we propose that the most reliable measurement crite- 
rion is the channel capacity itself. However, due to the fact 
that pursuing this basis requires excessive complexity in both 
calculations and time, we have just relied on the results of this 
section as a comparison measure throughout the paper 

Fig. [U illustrates the condition of 3 optimum matrices for an 
overloading factor of |, each optimized for a certain SNR, in 
an SNR sweep observation. As is clear, the matrices exhibit a 
channel capacity close to the extreme upper bound, i.e. 1, even 
in low SNRs. It is noteworthy to mention that the per-user 
channel capacity is defined as the ratio of the sum channel 
capacity to the number of users, keeping in mind that the 
channel is fairly divided among all users. 

The calculation cost of the capacity criterion is so immense 
that it cannot be extended to wide dimension ranges, thus the 
need for more simpler criteria comes to life. As previously 
explained, this is the main focus of this research. Henceforth, 
we have basically pursued a stagewise procedure in search 
of these alternate criteria, which immensely contribute to the 
possibility of designing highly reliable CDMA systems. 

In the next subsection, we have provided the Bit Error Rate 
criterion as a better-handled criterion than that of capacity. 
The BER, not only enables extensive search dimensions, but 
is observed to be reversely related to CDMA channel capacity. 

B. The BER Scenario 

In this subsection we have introduced another measure to 
find the optimum matrix, i.e. the Bit Error Rate. The optimality 
of a CDMA signature matrix could be based on both maximum 
channel capacity or minimum bit error rate. In our research, 
we are only concerned with the prior of the two standards, 
however, it seems that they both are closely bonded in a reverse 
relationship. By this, we conjecture that the lower the BER a 
certain signature matrix demonstrates, the higher its channel 
capacity will be. Using Genetic Algorithm to minimize the 
BER leads to an optimum matrix. Nevertheless, in order to 



implement a BER evaluation, we will have to utilize a channel 
decoder In fact, we have opted the ML0 decoder throughout 
this research. 

According to the aforementioned discussions, we carried out 
the same genetically based algorithm in order to determine 
optimum matrices with respect to minimum BER. 

It is apparent from Fig. [T] that our conjectures are in 
consent with reality. The minimum BER matrix adheres to the 
approximate upper bound for capacity, found in the previous 
section. This observation clarifies the fact that in order to 
be able to find higher-dimensioned code matrices, we could 
exploit the BER as a better measure than the exhaustive 
capacity computations. 

Although calculating the BER involves much less complex- 
ity than the capacity itself, it still lacks the feasibility to be 
implemented on practical orders of dimension. It is true that in 
order for the BER to be reversely correlated to the capacity, a 
large array of bits must be produced, coded, sent through the 
channel and decoded at the receiver The lower the number of 
evaluated bits, the more dissimilar BER and channel capacity 
tend to be. To be more precise, a single bit transmitted through 
an AWGN channel, needs to be compared with all the possible 
output constellation points to ensure ML detection. 

In the next section a faster and more practical set of 
standards have been presented, that slightly adjust the trade- 
off between optimality and implementability in favor of the 
latter 

C. Output Constellation Schemes 

The abovementioned calculations induces the motive to 
propose some less complex, but optimum criterion which 
could be easily implemented; hence we altered our viewpoint 
towards the constellation of noiseless output vectors in the 
image field. 

In this subsection, we focus on maximizing the capacity 
with respect to the constellation of noiseless output vectors. 
Throughout this paper we consider the output constellation 
points as, 

z.^AX, Xi,X2,--- ,X2. e {±1}"^! (10) 

At first, we focus on the minimum distance criterion, which 
is widely accepted among coding theory researchers. In this 
scheme, we have proposed that if we want to obtain an 
optimum matrix regarding to the capacity, we must maximize 
the minimum distance between the noiseless F's. Interestingly 
enough, the same measure, i.e. the minimum distance maxi- 
mization, yields sub-optimum results in our point of concern. 

The following relation explicitly explains this criterion, 

vi = min \\Zi — Zj\\ (11) 

where \\U\\ represents the Euclidean norm of the vector U. 
We shall denote vi as the Minimum Distance (MD) criterion 
throughout this text. Note that vi must be maximized in order 
to find a sub-optimum matrix. 

'^Maximum Likelihood 
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Fig. 2. Comparison Between the Sum Channel Capacity of the Different 
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Fig. 3. Comparison Between the Bit Error Rate of the Different Constellation 
Methods Matrices {ui, V2 and vz) and the Optimum BER Matrix Versus 
Different SNRs (/? = I) 



To take a more analytical approach, the upper bound for 
error probability is needed. Therefore the independency of x^'s 
results in. 



^e = ^E^( U \\Y.~Z,f>\%^Zu\A 
^ — 1 ^ I — 1 ' 



(12) 



k=l 



Since the noise vector elements are i.i.d Gaussian random 
variable with variance aN, the error probability upper bounded 

as. 



i=l 3 = 1 



\Zi - Zj 
2a N 



(13) 



where Q{x) is the Q-function, the Tail Probability of a zero 
mean Gaussian random variable with variance 1. 

According to (fTsT l. the Q Distance (QD) measure is intro- 
duced as follows. 
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(14) 



It must be considered that v^, as the QD criterion, should be 
minimized in order to find optimum code matrices. 

The 2-function that appears in V2 can be approximated with 
the less complex equation as below. 



Q(a;) 



aexp 




(15) 



where curve fitting results in a = 0.7, h—\ and c~ 1.6. 
According to (flST l, vi can be simplified as. 



\\z.-z,\ 



1 



1.6 



(16) 



where 1/3 refers to as the Exponential Distance (ED) criterion, 
which must be minimized. 



In high SNRs equation ( fT6b can be reduced to only one 
exponential element which corresponds to the Minimum Dis- 
tance; consequently we expect the similarity between MD and 
ED in high SNRs followed by the advantage of ED in low 
SNRs. 

In Figs. |2] and [3] the mentioned criteria, i.e. MD, QD, 
and ED, optimum matrices are compared together in separate 
capacity and BER plots. As can be seen, ED and QD have 
approximately the same performance, which demonstrates 
the perfect approximation in ( fTSl l. This provides the system 
designer with an advantage of less computational complexity 
in choosing perfect signature matrices alongside negligible 
behavioral offshoot. In addition, these two measures perform 
quite better than the simpler MD method, which can be sheerly 
proved due to the excessive complexity of the prior two 
schemes. 

In the next section, we have illustrated the optimality char- 
acteristics of the different proposed criteria, which is followed 
by in-depth comparisons between these schemes. The trade- 
off between matrix optimality and computational complexity 
is perfectly clear from analytical as well as simulation results. 

IV. Comparison and Results 

As an aftermath to the achieved criteria in the previous 
sections, we carried out their optimization via Genetic Algo- 
rithms. It is obvious that for each n, m and SNR, we expect 
a particular, and probably inconsistent, signature matrix to 
fulfill the optimality criterion. However, it was observed, and 
then certified, that optimum matrices in the same dimensions 
usually exhibit near-optimum, if not optimum, characteristics 
in other SNRs, hence providing perfect signature matrices for 
a wide extent of noise powers. Fig. |4] contains the per-user 
capacity versus SNR for the matrices optimized with capacity, 
BER, MD and ED criteria compared with a Welch Bound 
Equality (WBE) matrix set. As it is anticipated, because of 
optimality of WBE matrices for Gaussian-input distribution, 
there is no guarantee to be optimum as far as we are concerned. 
Note that the overloading limitation sets the extra restriction 
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Fig. 4. Comparison Between the Sum Channel Capacity of the Different 
Criteria Matrices Versus Different SNR's (/3 = f ) 

that such systems will fall short of the maximum feasible per- 
user capacity, that is 1, for low SNRs. 

It should be noted that in order to be able to compare the 
results of our criteria with the achievable capacity, we were 
obliged to evaluate the results in low orders of dimensions. 
However, due to the satisfying characteristics of the matrices 
found by optimizing these criteria, it can be predicted that 
optimizing these criteria in larger dimensions also provides 
near-optimum matrices. 

It is worth to mention that we have disregard the QD 
measure, due to exactly similar results to the simpler ED 
method. Considering Fig|4] it is clear that the optimum matrix 
for 8 dB is not optimum for other SNRs, especially high SNRs, 
so its falling below other matrices is not an issue. 

It can also be observed that the BER and ED criteria 
have relatively near optimum performance, even though the 
implementation is obtained with less complexity. Considering 
the fact that these two methods minimize the Bit Error Rate, 
i.e. error probability, the conjecture proposed in section |lll] 
is confirmed. Remarkably enough, despite the fact that ED 
and BER extremely differ in computational complexity, i.e. 
estimating the BER is more challenging due to the simula- 
tion of the channel and randomness behavior; ED and BER 
demonstrate the same performance regarding to sum channel 
capacity. 

MD, which is proved to be a specialization of the general 
ED, thusly falls beneath the other criteria in the figure. This 
sub-optimality may be thought to be compensated with the 
simplicity of MD. Of course, as anticipated, in high SNRs 
the MD and ED criteria tend to converge; providing an extra 
degree of freedom for the overloaded system. 

Another important feature of our proposed criteria-based 
technique is that it can lead to results with an extremely 
flexible behavior in different overloading factors (^), thus 
upheaving the implementability of overloaded CDMA sys- 
tems. Therefore, we have presented the results of the different 
criteria's per-user capacity while changing the overloading 
factor in Fig. |5] Such a feature brings about the extension 
of our approach's usage in limited available chip-rates, i.e. 
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Fig. 5. Per-user Capacity of Different Proposed Criteria Versus a Range of 
Overloading Factors. (In SNR 8dB) 

small m. It can be seen, the fact that we did not restrict our 
signature matrix elements in all the different scenarios, leads 
to us enabling high capacity quantities while the overloading 
factors are relatively high. Although the per-user capacity still 
decreases as the overloading factor grows, it descends more 
smoothly in our case; sustaining an acceptable range for larger 
overloads. 

V. Conclusion and Future Works 

In this paper, we proposed a method to determine the 
capacity related to a specific matrix, then in order to find 
optimum capacity matrices, a number of criteria have been 
proposed. These measures differ in optimality and complexity 
with one another. Generally, the proposed Exponent Distance 
(ED) criterion presents the best performance regarding the 
complexity and optimality. It should be stated that to find 
the optimum matrices in each criteria, we implemented a 
novel Genetic Algorithm approach. The real valued matrices 
found by this method discriminate our approach with the prior 
researches, providing outstanding capability in specific cases 
in which the chip rate is limited. Considering this feature, 
the optimum matrices tend to sustain an appropriate per-user 
capacities for high overloading factors. 

The algorithmic procedure in providing the optimum over- 
loaded signature matrices as well as precise analytical evalu- 
ations are appropriate future works in the aforementioned re- 
search area. Such algorithms may be utilized to gain even more 
time efficiency while designing an optimum overloaded code 
division multiple access system. In addition, implementing a 
similar approach for Optical CDMA systems in presences of 
Poisson noise is considerable. 
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